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Abstract 

A quantum kinetic equation is derived for the description of pair production in a time- 
dependent homogeneous electric field E(t). As a source term, the Schwinger mechanism for 
particle creation is incorporated. Possible particle production due to collisions and collisional 
damping are neglected. The main result is a kinetic equation of non-Markovian character. 
In the low density approximation, the source term is reduced to the leading part of the well 
known Schwinger formula for the probability of pair creation. We discuss the momentum and 
time dependence of the derived source term and compare with other approaches. 



1. Introduction 

Particle production in ultrarelativistic heavy-ion collisions raises a number of challenging 
problems. One of these inteijj|tjji||jj||}g|jons is how to incorporate the mechanism of-|j 



creation into a kinetic theoryH H U □ D Ullffi. Within the framework of a flux tube modelQLjOO, 
a lot of promising research has been carried out during the last years. In the scenario where a 
chromo-electric field is generated by a nucleus-^jpeleus collision, the production of parton pairs 
can be described by the Schwinger mechanisrrOuj'Ej. The produced charged particles generate 
a field which, in turn, modifies the initial electric field and may cause plasma oscillations. 
The interesting question i-af the back reaction of this field has been analyzed within a field 
theoretical approachli3'E£Hij. The results of a simple ^hcrhpm enorogical consideration based on 
kinetic equations and the field-theoretical treatmentuO'EjO agree with each other. The source 
terna-which occurs in such a modified Boltzmann equation was derived phenomenologically in 
RefE3. However, the systematic derivation of this source term in relativistic transport theory 
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is not yet fully carried out. Fof-ie»ample, recently it was pointed out that the source term may 
have non-Markovian characteioE3 even for the case of a constant electric field. 

In the present work, a kinetic equation is derived in a consistent field theoretical treatment 
for the time evolution of the pair creation in a time-dependent and spatially homogeneous 
electric field. This derivation is based on the Bogoliubov transformation for field operators 
between the asymptotic in-state and the uj simultaneous state. In contrast with phenomeno- 
logical approaches, but in agreement withElLJ, the source term of the particle production is 
of non-Markovian character. The kinetic equation derived reproduces the Schwinger result in 
the low density approximation in the weak field limit. 

The paper is organized as follows. In Section 2 we derive the kinetic equation with a new 
source term that describes particle production in a strong electric field within a field theoretical 
approach for both fermions and bosons. In Section 3 we discuss the main properties of the 
result and give numerical estimates for the case of a constant field and a model time dependent 
field. We summarize our results in Section 4. 

2. Dynamics of pair creation 



2.1. Creation of fermion pairs 



In this section we demonstrate the derivation of a kinetic equation with a source term 
for fermion-antifermion production. As an illustrative example we consider electron-positron 
creation, however the generalization to quark-antiquark pair creation in a chromoelectric field 
in Abelian approximation is straightforward. For the description of e + e~ production in an 
electric field we start from the QED lagrangian 



(1) 



where F^ v is the field strength, the metric is-jtaken as g^ u — diag(l, —1,-1,-1) and for the 7- 
matrices we use the conventional definition C3. In the following we consider the electromagnetic 
field as classical and quantize only the matter field. Then the Dirac equation reads 



(i7 M (9 M — e-y^A/j, — m)ip(x) = . 



(2) 



We use a simple field-theoretical model to treat charged fermions in an external electric field 
charactarized by the vector potential A^ — (0, 0, 0, A(t)) with A(t) = As(t). The electric field 

E{t)=E 3 {t) = -A{t) = -dA{t)/dt (3) 

is assumed to be time-dependent but homogeneous in space (E\ — E2 =0). This quasi- 
classical electric field interacts with a spinor field ifi of fermions. We look for solutions of the 
Dirac equation where eigenstates are represented in the form 



i-y°d + -y h p k - ey z A(t) + ■ 



(4) 



where k = 1,2,3 and the superscript (±) denotes eigenstates with the positive and negative 
frequencies. Herein R r (r = 1, 2) is an eigenvector of the matrix 7°7 3 



Ri 




R2 = 




(5) 



A quantum kinetic equation. . . 3 



so that R^Rs = 2<5 rs . The functions x (j>jt) are related to the oscillator-type equatl 



X (±) (p,i) = 



^(p,t)+ieA(t) ) x (±) (P,t) 



(6) 



where we define the total energy uj 2 (p,t) — e\ + Pa (t) , the transverse energy e\ = m 2 +p± and 
P\\(t) = p\\ — eA(t). The solutions X^iPit) of Eq. (g) for positive and negative frequencies 
are defined by their asymptotic behavior at to — t — > — oo, where A(to) — 0. We obtain 



X (±) (P, *) ~ exp ( ± iwo (p) t) , 



(7) 



where the total energy in asymptotic limit is given as uio(p) — uo(p,to) = lim uj(p,t). 

l—l * — > — oo 

Note that the system of the spinor functions mft is complete and orthonormalized. The field 
operators ip(x) and 4>(x) can be decomposed in the spinor functions m) as follows: 



r,p 



^ pr \x) bpr(to) + ^t\x) dt pr {to) 



(8) 



The operators bp r (to), b pr (to) and d pr {to) , d pr (to) describe the creation and annihilation of 
electrons i-apd positrons in the in-state |0i n > at t = to, and satisfy the anti-commutation 
relations Eil 

{MfcO.&ivfa)} = {dpr(t ),d+ r ,(t )} = S rr , 8 pp , . (9) 
The evolution affects the vacuum state and mixes states with positive and negative energies 
resulting in non-diagonal terms that are responsible for pair creation. The diagonalization of 
the hamiltonian corresponding to a Dirac-particle (Eq. (□)) in the homogeneous electric field 
(^) is achieved by a time-dependent Bogoliubov transformation 



b pr (t) = a p (t) bp r (to) + pp(t) d± Pr {to) , 
d pr {t) = a-p(t) d P r(t ) - P-p{t) bt pr (t ) 



with the condition 



m*)i 2 + mt)\ 2 = i 



(10) 



(11) 



Here, the operators b pr (t) and d pr (t) describe the creation and annihilation of quasipar- 
ticles at the time t with the instantaneous vacuum 0t >• Clearly, the operator system 
b(to), b + (to); d(to), d + (to) is unitary non-equivalent to the system b(t),b + (t);d(t),d + (t). The 
substitution of Eqs. ([jlj into Eq. (^) leads to the new representation of the field operators 



ip{x) 



E 

r,p 



*£\x) M*) + dt pr (t) 



(+)/ 



(12) 



The link between the new 9^P(x) and the former (^) basis functions is given by a canonical 
transformation 



4~\x) = a p (t) ¥ p -\x)-l3 P (t) *<+>(*) , 
^\x) = a P (t) yM(x) pr + 0p(t) *£\x) 



(13) 
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Therefore it is justified to assume that the functions have a spin structure similar to that 

in Er, ("Eh 



of in Eq. 



^pr (-^) 



17° do + "f k pk — e7 3 A(t) + m 
where the dynamical phase is defined as 

®{p,t) = / dt'u)(p,t') 



i ±] (x) R r e 



ii©(£) ^ipx 



(14) 



(15) 



The functions <j>^ are yet unknown. The substitution of Eq. (^) into Eqs. (^) leads to the 
relations 



X (+) (p,i) = «;(t)^(t)e 



+ 4"»(t) 



(16) 



-i&(P,t) 



Now we are able to find explicit expressions for the coefficients ap(t) and /3p(t). Taking into 
account that the functions X^(Pi t) are defined by Eq. (^|), we introduce additional conditions 
to Eqs. (|l^) according to the Lagrange method 



-iw(p, t) 



a p (t) 4><f\t) e-' 6 ^" + /31(f) 4 + \t) e ie <^> 



(17) 



X {+) {P,t) = iu{p,t) 



Mm P i&(P,t) 



«;-(*) ^y\t) e 



Pp(t) <f> ( i\t) 



-i&(p,t) 



Differentiating these equations with respect to time, using Eqs. (g) and (Jl|) and then choosing 
as an Ansatz 



we obtain the following differential equations for the coefficients 



(18) 



<%(f) = 



eE(t)e± a * uy aie(p,t) 



2^ 2 (p,t) 

eg(t)ex 
'2io 2 (p,t) 



ap(t) e- 2 ^^ 



(19) 



As the result of the Bogoliubov transformation we obtained the new coefficients of the instan- 
taneous state at the time t. The relations between them read 



>*(«) = 



/?(*) = 



2^0), t) (w(p,f)-P||(f)) 



2y/u(p,t) (w(p,t)-Pj|(t)) 



fw(p,t) X H fet)-iX H fet)) e" ie ^ 



(20) 
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It is convenient to introduce new operators which absorb the dynamical phase 

B pr (t) = b pr (t) e - ie <^) , (21) 
D pr (t) = d pr (t) e" 1 ®^*) (22) 

satisfying the anti-commutation relations: 

{B pr (t),B+ r ,(t)} = {D pr (t),D+ r ,(t)} = S rr , S pp , . (23) 

It is easy to show that these operators satisfy the Heisenberg-like equations of motion 

dB pr (t) _ eE(t)e±_ 



dt 2co 2 (p,t) - p 

dD pr (t) eE{t)e 



D± pr (t)+i [H(t), B pr (t)] , 
Bt p M+i Wit), D pr (t)] , 



(24) 



dt 2uj 2 (p,t) ~~ v 

where H(t) is the hamiltonian of the quasiparticle system 



H(t) =Y,u>(p,t)l B+ r {t) B pr {t) - D. pr (t) D+ pr (t) . (25) 



The first term on the r.h.s. of Eqs. (|24|) is caused by the unitary non-equivalence of the 
in-representation and the quasiparticle one. 

Now we explore the evolution of the distribution function of electrons with the momentum 
p and spin r defined as 

fr(p,t) =< in |6^.(i) b pr (t)\O tn >=< Oi„\B±(t) B pr {t)\O ln > . (26) 

According to the charge conservation the distribution functions for electrons and positrons are 
equal f r (p,t) = f r {p,t), where 

f r (jp,t) =< O in \d+ pr (t) d-p r (t)\O in > = < in \ D+ Pr (t) D- Pr (t)\O in > . (27) 

The distribution functions ( pi|) and (^) are normalized to the total number of pairs N(t) of 
the system at a given time t 

J2 MP, t)=J2 f'fc *) = N ® ■ ( 2§ ) 
r,p r,p 

Time differentiation of Eq. ( ptj ) leads to the following equation 

df r (p,t) eE{t)e A 



dt u 2 (p,t) 



Re{$ r (p,t)} . (29) 

Herein, we have used the equation of motion (|24j) and evaluated the occurring commutator. 
The function & r (p, t) in Eq. (^) describes the creation and annihilation of an electron-positron 
pair in the external electric field E(t) and is given as 

$r(p,t) =< O ln \D- pr (t) B pr {t)\O m > . (30) 
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It is straightforward to evaluate the derivative of this function. Applying the equations of 
motion we obtain 



d$ r (p,t) _ eE(t)s s 



dt 



2tu 2 (p,t) 



2/r (P,t) 



2iu{p,t) $ r (p,t) 



(31) 



where because of charge neutrality of the system, the relation f r (p,t) = f r (p,t) is used. The 
solution of Eq. (^l|) may be written in the following integral form 



ex 



dt 



, eE(t') 
w»(p,tf) 



2i[B(p,t')-B(p,t)] 



(32) 



The functions 0(p, t) and Q(p,t') in Eq. (|32j) can be taken at to (see the definition (|15[)). 
Hence with A(to) = 0, the function $ r (p, t) vanishes at to- Inserting Eq. ( p2[ ) into the r.h.s of 
Eq. (S) we obtain 



df r {p,t) _ eE(t)e A 



dt 



2w 2 (p,t) 



,eE{t')e A 



2/ r (p,t') 



oos 2[e(p-,t)-0(p,t')] 



(33) 



Since the distribution function obviously does not depend on spin (^3|), we can define: f r = 
/. With the substitution f(p,t) — > F(P,t), where the 3-momentum is now defined as 
P(p±, P\\ (t)), the kinetic equation (p3[) is reduced to its final form: 



dt 

with the Schwinger source term 



dF(P,t)_8F(P,t) anP,ii =Si 



>(-) 



eE(t)ex 
2co 2 (p,t) 

A 



dt 



dt 



,eE{t')e A 



l (p,t') 



dPn (t) 



1 - 2F(PA') 



(-) 



cos [2[e(p,t)-e(p,t')] 



(34) 



(35) 



Recently in RefE, a kinetic equation similar to ( p4| ) has been derived within a projection 
operator formalism for the case of a time-irzdependent electric field where it was first noted 
that this source term has non-Markovian character. As well as in this method the multiple 
pair creation is not considered. The presence ofptiie Pauli blocking factor [1 — 2F(P,t)] in 
the source term has been obtained earlier in RefEfl. We would like to emphasize the closed 
form of the kinetic equation in the present work where the source term does not include the 
anomalous distribution functions (J3C|) for fermion-antifermion pair creation (annihilation). 



2.2. Creation of boson pairs 

In this subsection we derive the kinetic equation with source term for the case of pairs of 
charged bosons in a strong electric field. 
The Klein-Gordon equation reads 



(9 M + ieA>*)[da + ieA^) + ■ 



(x) = 0. 



(36) 



The solution of the Klein-Gordon equation in the presence of the electric field defined by the 
vector potential = (0, 0, 0, A(t)) is taken in the form Ej 



^\x) = [2w(p)] 



-i/2 jap (±) 



(P.*) 



(37) 
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where the functions g ( - ± \p,t) satisfy the oscillator-type equation with a variable frequency 

g {± \p,t)+uj 2 {p,t) 5 (±) (p,t)=0 . (38) 

Solutions of Eq. for positive and negative frequencies are defined by their asymptotic 

behavior at to = t — > — oo similarly to Eq. (Q) . 
The field operator in the in-state is defined as 

<P(x) = Jd 3 p[ 0^ (x) a p (t ) + 4 +) (x) 6+ (to) ] . (39) 

The diagonalization of the hamiltonian corresponding to the instantaneous stateis achieved by 
the transition to the quasiparticle representation. The Bogoliubov transformation for creation 
and annihilation operators of quasiparticles has the following form 

a p (t) = ap(t) Op (to) + (3p(t) b± f (to) , 

(40) 

b- p (t) = a- p (t) bp(t ) + P- P (t) a± P (t ) 

with the condition 

K(*)| 2 - Wp(t)\ 2 = 1 • (41) 
The derivation of the Bogoliubov coefficients a and (3 is similar to that given in the previous 
subsection. We obtain the equations of motion for the coefficients of the canonical transfor- 
mation ( [io| ) as follows 

Mt) - (42) 

0p{t) = M^T) a f {t)e ■ (43) 

Following the derivation for the case of fermion production, we arrive at the final result 
for the source term in the bosonic case 



p _ eE(t)s ± [' , eE{t')e ± 



l + 2F(P,t') 



cos [2[e(p,t)-Q(p,t')] ), (44) 



which differs from the fermion case just by the sign in front of the distribution function due 
to the different statistics of the produced particles. 

3. Discussion of the source term 



3.1. Properties of the source term 

We can combine the results for fermions (|3f]) and for bosons (^) into a single kinetic 
equation 

dF, ±) (P,t) dF (±) (P,t) , BF (±) (P,t) . — . 

Here, the upper (lower) sign corresponds to the Bose-Einstein (Fermi-Dirac) statistics. Based 
on microscopic dynamics, these kinetic equations are exact within the approximation of a 
time-dependent homogeneous electric field and the neglect of collisions. The Schwinger source 
terms ( |3q ) and (jlij) are characterized by the following features: 
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1. The kinetic equations ( |45| ) are of non-Markovian type due to the explicit dependence 
of the source terms on the whole pre-history via the statistical factor 1 ± 2F(P, t) for 
fermions or bosons, respectively. The memory effect is expected to lead to a modification 
of particle pair creation as compared to the (Markovian) low-density limit, where the 
statistical factor is absent. 

2. The difference of the dynamical phases, Q(p,t) — <d(p,t'), under the integrals ( pS5| ) and 
([lij) generates high frequency oscillations. 

3. The appearance of such a source term leads to entropy production due to pair creation 
(see also Raua) and therefore the time reversal symmetry should be violated, but it does 
not result in any monotonic entropy increase (in absence of collisions). 

4. The source term and the distribution functions have a non-trivial momentum dependence 
resulting in the feat that particles are produced not only at rest as assumed in previous 
studies, e.g. RefJiSI. 

5. In the low density limit and in the simple case of a constant electric field we reproduce 
the pair production rate given by Schwinger's formula 

scl =zrJ^ r39 J d3ps ^ t)=e ^^{-w\) • (46) 

As noted above, Raul found that the source term has a non-Markovian behaviour by 
deriving the productipin rate within a projector method. In the limit of a constant field our 
results agree withRefR In our approach the electric field is treated as a general time dependent 
field and hence there is no a priori limitation to constant fields. However our result allows 
to explore the influence of any time-dependent electric field on the pair creation process. It 
is important to note that in general this time dependence should be given by a selfconsistent 
solution of the coupled field equations, namely the Dirac (Klein-Gordon) equation and the 
Maxwell equation. This would incorporate back reactions as mentioned in the introduction. 
However, the solution of such a system of equations is beyond the scope of this work. Herein 
we will restrict ourselves to the study of some features of the new source term. 

Finally we remark that the source term is characterized by two time scales: the memory 
time 

e_L 

Tmem ^ \4( ) 

eE 

and the production interval 

Tprod = V < S(±) > , (48) 

with < > denoting the time averaged production rate. As long as E << m 2 /e < e\/e, 
the particle creation precess is Markovian: T mem « r vro d- This results for constant fields 
agree with those of Rautj. 



3.2. Numerical results 

In order to study the new source term, we consider two different cases, namely a constant 
field and a time dependent electric field. For the numerical evaluation we start with Eq. ( pB| ) 
assuming a dilute system, F = 0, and 
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30.0 




-3.0 -1.0 1.0 3.0 

P,/e x 

Fig. 1. The pair production rate as a function of parallel momentum for a constant, strong field 
(upper plot: Eq = 1.5) and a weak field (lower plot: Eo = 0.5) at t = 0. 



where we have introduced dimensionless variables 

E(t) = eE(t)/e 2 ± , i = tej_ , (50) 

P\\ =P||/e_L, U) = u/e±. (51) 

This notation is convenient to distinguish the weak field (E < 1) and strong field [E > 1) 
limits. A particularly simple result is obtained if we assume a constant field, 

A(i) = A(i)/e± = tEo/e , (52) 

where the electric field does not depend on time and the energy is given as 

w£(P||,*) = l + (P|| -E t) 2 . (53) 

For the source term we obtain 

5(P|, , t) = fl I dt' - ! cos ( 2 f dt"ti (P\l , t")) ■ (54) 

In Fig. ^ we plot the particle production rate as function of the parallel momentum py for 
a weak constant field and a strong field, respectively. The rates are normalized to be of the 
order of one for Eo — 0.5 at zero values of both momentum and time. The production rate 
is positive when particles are produced with positive momenta. Pairs with negative momenta 
are moving against the field and hence get annihilated, clew-ly to be seen in the negative 
production rate. These results mainly agree with those of Ref.Q, since no time dependence was 
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■10.0 

1.0 

_i 0.0 




-3.0 2.0 7.0 

K 

Fig. 2. The pair production rate as a function of time for a constant, strong field (upper plot: 
Eg = 1.5) and a weak field (lower plot: Eq = 0.5) at zero parallel momentum. 



considered. Note that using thft prefactor and field strengths chosen by Rau, we reproduce 
exactly the results given in Ref.El. 

6 considering the time dependence of the source term, we go beyond the analysis of 
In Fig. we display the time dependence of the production rate at zero momentum. 
The maximum of the production rate is concentrated around zero and shows an oscillating 
behaviour for large times. Indeed it is possible to write Eq. (JH^) in terms of the Airy function 
because the constant field provides an analytical solution for the dynamical phase difference 
using the energy given in Eq. (|H^). The production of pairs for strong fields is larger than 
that of weak fields, and the typical time period of the oscillations is smaller. 

The situation changes if we allow the electric field to be time dependent. We assume a 
Gaussian field at the dimensionsless time r with the width a = ae±, 



E(t) = E e 



-(i-r) 2 /5 



(55) 



and obtain 



A(t) = -E a- 



Erf[(t - T)/a) - Erf[(-T - t ) /a] 



(56) 



The occuring error function is defined as 



Erf(z) 



x dx . 



(57) 



Using this Ansatz for the field strength in Eq. (jf9|) we obtain the numerical results plotted 
in Figs. ^ and ^. Therein all occuring values are dimensionsless. The electric field is non 
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Fig. 4. The pair production rate, S(p|| , t), as a function of time and parallel momentum for a time 
dependent strong electric field charactarized by Eq = 1.5, a = 1 and r = . All plotted values 
are dimcnsionless as described in the text. 
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zero within a certain width a = 1 for r — around t = 0. Therefore the oscillations for 
times beyond the time where the electrical field is finite are damped out. The pair production 
rate is peaked around small momenta for both weak and strong fields. For strong fields the 
distribution is shifted to positive momenta but remains still close to small parallel momenta. 
It is important to point out that the production of particles happens not only at,rest (py = 0) 
what is assumed in many works also addressing the back reaction problem, e.g.Ej. In contrary 
we find a non-trivial momentum dependence of the pair creation rate depending on the field 
strength and on time. 

4. Summary 

We have derived a quantum kinetic equation within a consistent field theoretical treatment 
which contains the creation of particle-antiparticle pairs in a time-dependent homogeneous 
electric field. For both fermions and bosons we obtain a source term providing a kinetic equa- 
tion of non-Markovian character. The source term is characterized by a pair production rate 
that contains a time integration over the evolution of the distribution function and therefore 
involves memory effects. 

For the simple case of a constant electric field in low density limit and Markcwian approx- 
imation, we analytically and numerically reproduce the results of earlier worksu. Since our 
approach is not restricted to constant fields, we have explored the dependence of the source 
term on a time dependent (model) electric field with a Gaussian shape. Within these two 
different Ansatze for the field, we have performed investigations of the time structure of the 
source term. The production of pairs does not happen at rest only. We observe a non-trivial 
momentum dependence of the source term depending on the field strength and on time. 

The particle production source is dominated by two time scales: the memory time and the 
production time. The numerical results mainly show oscillations due to the dynamical phase 
and urge the need to include the Maxwell equation to determine the electric field by physical 
boundary conditions (back reactions), work in this direction is in progress. Furthermore, it 
would be of great interest to extend this approach to the QCD case to explore the impact 
of a non-Markovian source term on the pre-equilibrium physics in ultrarelativistic heavy-ion 
collisions. 
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